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Low loss Bulk Acoustic Wave devices are considered from the point of view of the solid state approach as 
phonon-confining cavities. We demonstrate effective design of such acoustic cavities with phonon-trapping 
techniques exhibiting extremely high quality factors for trapped longitudinally-polarized phonons of various 
wavelengths. Quality factors of observed modes exceed 1 billion, with a maximum Q-factor of 8 billion and Q X 
f product of 1.6 • 10 18 at liquid helium temperatures. Such high sensitivities allow analysis of intrinsic material 
losses in resonant phonon systems. Various mechanisms of phonon losses are discussed and estimated. 

AFabry-Perot interferometer is an optical resonant system that typically consists of two parallel highly 
reflecting mirrors. Due to their high quality factors, a Fabry-Perot cavity is often used to store photons for 
almost milliseconds while they bounce between the mirrors. Low loss is necessary in order to increase 
interaction between light and matter or other physical substances under study. Such systems are widely used in 
telecommunication, lasers, gravitational wave detection, spectroscopy, etc. 

Similar principles can be applied to the field of acoustics, bulk phonons may be resonantly confined in a thin 
film or a crystal lattice to build a phonon cavity. In this case, resonant phonons are stored between borders that 
separate different media playing a role of mirrors. 

Losses in acoustic resonators can be classified into two groups: first, losses attributed to intrinsic material 
processes, second, losses due to resonator design. The latter are mostly due to energy leakage through clamping, 
substrates, electrodes, etc 1 . The goal of the cavity design is to make "engineering losses" negligible in order to 
achieve a regime where losses are due to the material itself. If this regime is achieved the cavity quality factor can be 
increased by orders of magnitude by cooling to cryogenic temperatures. 

The thin film approach (realised in Thin Film Bulk Acoustic Resonators and High Overtone Bulk Acoustic 
Resonators) provides technology that can guarantee smaller distances between the borders, and thus higher 
resonant frequencies. The main disadvantages of these systems are the need for a substrate and the inability to 
efficiently control the film geometry due to limitations of the deposition process. Both result in additional losses 
due to phonon leakage into the environment. 

The crystalline lattice approach utilises structures built from pieces of single crystal material (e.g. quartz Bulk 
Acoustic Wave (BAW) resonator). Due to a completely different fabrication process, this technology can provide 
a more sophisticated design that guarantees sufficient limitation of phonon leakage. In addition, extremely low 
losses of such devices extend the possible frequency range to comparatively high frequencies. In this work, we 
show that this approach allows Q-factors of above 1 billion with a maximum of 8 billion corresponding to a Q Xf 
product of 1.6 • 10 18 at liquid helium temperatures. This is 310 times higher than the previous value achieved at 
this temperature 2 . Such high sensitivity should allow future applications of acoustical cavities to be broadened and 
be used as an alternative (to optical cavities) experimental base for broad range of physical experiments, including 
material characterisation 3 , quantum cavity electromechanics 4 ' 5 , quantum metrology 6 ' 7 etc. 

Results 

Phonon excitation and trapping. Vibrations of ideal plate resonators have been studied for many decades 8 . 
Unlike their idealised counterparts, actual cavities pose two almost contrary requirements. First, they have to be 
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Plate Disk Phonon Distribution 



Figure 1 | Cross section in the r — z plane of the quartz Bulk Acoustic 
Wave (BAW) resonator of cylindrical geometry with typical distribution 
of trapped phonons in the radial direction. 

coupled to the environment for effective excitation and measure- 
ments, second, they must be uncoupled from the environment in 
order to avoid phonon leakage. The solution to this corundrum is to 
excite vibrations with the electro -magnetic field while mechanically 
separating the crystal from the environment. 

Coupling to electromagnetic measuring and excitation fields is 
introduced by making a vibrating body from piezoelectric material. 
Imposing a potential difference between two surfaces of the plate 
creates a gradient of the electric potential along the plate thickness. 
The spatial Fourier harmonics of this gradient are coupled to mech- 
anical vibrations through the piezoelectric effect. Frequencies of 
these vibrations are set by plate dimensions, material properties 
and the order of the harmonic. This type of vibration is referred as 
a thickness vibration that has to be distinguish from parasitic flexure 
vibration modes. Detection of the plate vibration is performed by 
measuring the displacement current between the surfaces of the 
plate. 

Due to the surface boundary conditions, such a cavity can sustain 
resonant phonons of discrete wavelengths. When describing an 
acoustic mode in a piezoelectric plate, the number of half- waves in 
the plate thickness is referred as the overtone (OT) number, and the 
resonance itself as an overtone. The cavity geometry allows only anti- 
symmetric types of vibration to be excited piezoelectrically, thus 
limiting experimentally observed OTs only to odd numbers. 

In order to avoid any energy leakage, the vibrating part of the plate 
has to be mechanically separated from the environment. This is 
achieved by two means: first, the excitation electrodes are separated 
from the crystal plate by a tiny vacuum gap; second, acoustic energy 
is trapped in the centre of the disk plate while the disk is clamped 
from the sides. The trapping is done by introducing curvature into 
one of the plate surfaces by manufacturing a plano-convex shape 910 . 
A schematic cross-section of the cavity is shown in Fig. 1. This 
variable plate thickness increases the local phonon frequency along 
the plate disk radius closer to the edges. The corresponding energy 
difference along the radial direction due to this curvature can be 
considered as an effective potential well. So, phonons are trapped if 
their total transverse energy (corresponding to in-plane wave num- 
bers) do not exceed the energy difference between the edge and the 
centre. This situation significantly decreases the probability to find a 
resonant phonon near the clamping points. Theoretical treatment of 
the energy trapping problem from the point of view of linear elasto- 
dynamics is given in 9 . The calculations predict the mode shape in the 
form of polynomial-normal distribution. This result has been con- 
firmed by means of X-ray spectroscopy 11,12 . 

Details of cavity design. Tested acoustic cavities are made from Top 
High Quality (THQ) alpha synthetic quartz produced by Gemma 
Quartz & Crystal according to a hydrothermal growth process. 
Impurities are removed from the heart of the quartz bars before 
the latter are cut up in oriented blanks. Disks come from slices cut 
out from a synthetic quartz crystal bar according to specific angles 
with respect to the crystalline axis (so-called singly rotated cuts or 
doubly rotated cuts). The tested cavities are doubly- rotated Stress 




Figure 2 | Photograph of the electrode and cavity mounting arrangement 
of the quartz BAW cavities. 



Compensated (SC) cut plates 13 , which is one of the most commonly 
used cuts. Doubly rotated cut cavities support three possible phonon 
polarizations with different velocities: the quasi-longitudinal (V L = 
6757 m/s), fast quasi-shear (V FS = 3966 m/s) and slow quasi-shear 
(Vss — 361 1 m/s). All three types of polarisation are not pure due to 
crystal anisotropy. According to the contoured-resonator axis [x, y, 
z] 14 , the y-axis being normal to the plate, their polarizations are 
respectively: [0.226, 0.968, 0.111], [-0.198, -0.0657, 0.979], 
[0.957, -0.243,0.177]. 

After the slices are cut they are then rounded, lapped and polished 
resulting in plano-convex contoured plates. Polishing is a key step of 
the process: it consists of a mixing of mechanical polishing and 
chemical etching. This last operation is mandatory to remove the 
"dead layer", i.e. a damaged layer - a blend of amorphous and poly- 
crystalline randomly-oriented matter resulting from the first process. 
Chemical etching is also one efficient way to control precisely the 
thickness of the resonator, a key parameter for defining the frequency 
of waves propagating along the thickness axis. At least, mechanical 
strains induced by machining should also be removed by an addi- 
tional reheat. 

The cavity design for which both electrodes are not deposited on 
the vibrating plate itself is known as (BVA) 15 . Electrodes are depos- 
ited on two quartz supports placed on each side of the resonator, and 
in such a way that the gap between electrodes and resonator surfaces 
is typically between 5 and 10 ^m 210 . A photograph of an assembled 
cavity is sown in Fig. 2. 

In this work we compare performance of the following types of 
acoustical cavities: 

1. 1.08 mm thick, 13 mm diameter electrode-separated disk cav- 
ities initially designed to sustain shear vibration of 5 MHz at 
room temperature (manufactured by BVA Industrie) 2 ; 

2. 0.54 mm thick, 13 mm diameter electrode-separated disk cav- 
ities initially designed to sustain shear vibration of 10 MHz at 
room temperature (manufactured by BVA Industrie) 10 ; 

3. 1 mm thick, 30 mm diameter electrode -separated disk cavities 
with higher grade surface polishing initially designed to sustain 
shear vibration of 5 MHz at room temperature (manufactured 
by Oscilloquartz SA); 

4. 0.54 mm thick, 13 mm diameter electrode-deposited disk cav- 
ities initially designed to sustain shear vibration of 10 MHz at 
room temperature (manufactured at FEMTO-ST). 

It has to be noted that all cavities are initially designed for room 
temperature operation with slow- speed thickness -shear vibrations 
and not for longitudinally-polarized ones. This is due to the fact that 
the frequency-temperature characteristic of the former exhibits a 
turnover point, which is not the case of the latter. Moreover at liquid 
helium temperatures, longitudinal phonons exhibit significantly 
lower losses than any shear phonons 2 . Direct cavity optimisation 
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for liquid helium operation is still impossible due to unknown tem- 
perature coefficients of the elastic constants at these temperatures. 

It has to be noted that two parallel excitation electrodes form a 
parasitic capacitor (typically close to 3 pF) creating a parallel photon 
branch. At higher frequencies, the number of these non-resonant 
photons exceeds the number of acoustic phonons decreasing the 
mode contrast. Hence, the parasitic capacitance uncouples the acous- 
tic motion from the electromagnetic environment. Nevertheless, the 
excitation of phonons of shorter wavelengths is possible with effec- 
tive electrical compensation of this shunt capacitance. 



Discussion 

The four cavities described above were characterised for resonant 

longitudinal phonons of half wavelengths varying from - to — of the 

plate thickness. Q-factors were extracted by fitting measured electric 
admittance of the resonator around each resonance frequency. Fig. 3 
shows an example of experimental results for the 3rd type of acoustic 
cavities measured for the 5th OT at 3.75 K. 

The results (see Fig. 4) show that the cavity with the largest dia- 
meter and better polishing exhibits significantly lower losses. This 
fact could be attributed to better trapping (due to larger potential 
well) and lower scattering on surface imperfections. In addition, all 
four samples exhibited an increase of Q as a function of OT number, 
at least for a certain range of frequencies. This is explained by 
improving of the phonon trapping with increasing OT number as 
follows from the theory 9 . Temperature dependence of the selected 
overtones for this cavity is shown in Fig. 5. 

Total intrinsic losses of an acoustic cavity apart from lack of trap- 
ping are a sum of losses due to different mechanisms: 
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where the Q-values represent different mechanisms. These mechan- 
isms are explained below combining the results of frequency (Fig. 4) 
and temperature (Fig. 5) dependences. 
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Figure 4 | Measured quality factors for different OTs of the (quasi-) 
longitudinal vibration of the four cavities under investigation at 3.75 K. 
Error bars for standard deviations lower than ±2.5% are not represented. 

Both higher values in brackets are close to the measurement set-up 
limitation. 
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Figure 5 | Temperature dependance of cavity losses for different OTs of 
the (quasi-)longitudinal vibration in the cavity number 3. 

First, the most important mechanism of energy loss in acoustic 
cavities at low temperatures and relatively high frequencies is loss 
mechanism due to phonon-phonon scattering over crystal lattice 
anharmonicity (Qphonon-phonon in equation (1)). This mechanism 
was first analytically considered by Landau for the case of shear- 
polarized phonons 16,17 . It was also pointed out that longitudinal pho- 
nons do not exhibit this type of losses in the case of two phonon 
collision, due to significantly higher sound velocity. Later the theory 
was extended to the case of three colliding phonons 18 predicting the 



following law: Q 



r -l 

phonon-phonon 



T . This law is experimentally 



observed for the large acoustic cavity between between 6 K and 
13 K (Fig. 5). According to the theory 18 , the acoustic wave attenu- 
ation a is proportional to the resonance angular frequency cd as long 

as -^-^ < < 1, and cot^M where t denotes the average lifetime of 

thermal phonons, h the Planck constant, and k B the Boltzmann con- 
stant. As a consequence, the Qphonon-phonon-fa&ox °f an acoustic wave 
propagating with velocity V 

nw , N 

Q= — 2 

would no longer depend on frequency co. Since these conditions are 
met in the present experiment, the Q-factor tends to be independent 
off(Q = const law) as shown in Fig. 4. 

Second, for temperatures below 6 K the Q(T) dependence changes 
the slope tending to T"° 3 scaling law (Fig. 5). Such dependence have 
been observed in a various of electromechanical systems 19 " 23 , and 
acoustic cavities at lower temperatures 10 , in particular. The T" 1/3 
dependence is theoretically explained by coupling to so-called "two 
level systems" (TLS) 24 (Qtls in equation (1)). In synthetic quartz 



these TLSs are attributed to impurity ions, typically A/ 3+ , Na + , Li + , 
Si 4+ and H + , etc 25 " 28 studied by electron-spin resonance, infrared 
spectroscopy, dielectric relaxation spectroscopy and thermolumines- 
cence. Although the utilised material is premium pure, concentration 
of these ions could reach 1-3 ppb. 

Third, the increase of cavity losses at higher frequencies occurs for 
the most of resonators (see Fig. 4), and can be explained by phonon 
scattering on the cavity surfaces {Qscattering in equation (1)). 
Scattering originates from the surface roughness as well as from 
the presence of a dead layer. Indeed, the reflection coefficient on a 
rough surface with a roughness variance cr 2 , assuming a gaussian 
probability density, and can be expressed as 29 R = e~ 2k2(j2 , where k 
is the wave number. At the same time, propagating through the 
thickness /z, the acoustic wave is attenuated by e~ xh . Thus, the 
reflection is equivalent to a modified path causing an attenuation 
equal to a' = a + ln(R/h). Consequently, for an OT number n such 
that kh = nn, the associated Q-factor becomes: 



h 2 



{.scattering z 



2n<i 2 



(3) 



Since the resonant frequency w is almost proportional to the over- 
tone number n, this mechanism results in Q X w = const law. The 
surface scattering mechanism strongly depends on the polishing 
quality denoted by the roughness variance a. This explains the better 
performance of cavity number 3 (Fig. 5), which has a better surface 
polishing grading. The largest observed values of quality factor were 
at n = 65 for cavity number 3, which has a a of one nm (consistent 
with specification of Oscilloquartz fabrication process). For this type 

2 

of cavity, Q sca ttering is estimated to be - • 10 , which gives an upper 

limit of 3 X 10 10 for the highest measured OT n = 65 at 204 MHz. 

For other cavities the polishing process is inferior with a varying 

from 4 to 7 nm. For example, for cavity type 2 with a specified a 

g 

= 4 nm, Qscattering is found to be - • 10 9 and gives a value of 4.2 X 10 8 

atrc = 19(118. 6 MHz), which is consistent with the experimentally 
measured value of 3.3 X 10 8 . Thus, we can conclude that the highest 
measured Q-values are close to estimated values of Qscattering* dem- 
onstrating that the cavity roughness could be a major factor limiting 
the Q-factor for high frequencies. 

Fourth, although thermoelastic losses do not exist for pure shear 
vibration, since there is no volume change, they are present for the 
longitudinal modes. So, the Qthermoeiastic for longitudinal vibration 
can be estimated based on the description of the wave propagation 
including the coupled equations of motion of the medium and ther- 
mal conduction 30 : 



PUi = Cij k iUk,lj- 



-pc s r- 



T 0 XkiUk,i = 0 



(4) 
(5) 



with Xij = Cfjuotkb an d where u { is the displacement vector, c ijk i the 
tensor of isothermal elastic constants, a^- the thermal expansion 
tensor, the thermal conductivity tensor, p the volumic mass, 
and C s the specific heat at constant strain, 3T the temperature 
change around the reference temperature T 0 . Although, most of 
these coefficients for the lightly anisotropic quartz are unknown at 
cryogenic conditions, simplified calculations can be made assum- 
ing an isotropic medium and propagation along the Z-axis. The 
optical Z-axis is the most experimentally studied case. This ana- 
lysis implies calculation of a turning-point frequency /* that 
shows the regime of acoustic wave propagation - adiabatic or 
isothermal: 



r 



27% 



(6) 
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where 



X3- 



for 



— ^ is the material diffusivity, and V 3 = , 
pC s y P 

propagation along the Z-axis. The former is calculated based on 
data on thermal conductivity and heat capacity at cryogenic tem- 
peratures 31 ' 32 giving 1 < /* < 5.5 MHz for 3 < T < 10 K. Over 
this range, all the resonant frequencies are always greater than/*, 
and acoustic wave propagates in the quasi-isothermal regime. For 
this regime, the maximum attenuation is: 

_ e 3 nf* 



Ql 



where Q L is a mode loaded quality 



(7) 



with 63 -- 



, where X 33 = 2c u (x u + c 33 ot 33 , because a lx = a 2 2 



pC% 3 

and c 13 = c 23 in case of crystalline quartz. Whereas the thermal 
expansion coefficients have been measured at low temperatures 31 , 
the elastic constants on the other hand can just be roughly 
approximated by their values at room temperature. Nevertheless, 
the corresponding thermoelastic losses estimated from equation 
(2) remain always at least two orders of magnitude lower than 
those measured below 10 K. Since aoo is almost independent of 
frequency for a given material, these type of losses result in Q/co = 
const law. 

We demonstrate that acoustic cavities with effective trapping of 
longitudinally polarised phonons can exhibit quality factors above 1 
Billion at liquid helium temperatures and frequencies 10-200 MHz. 
This shows that phonon cavities can be considered as an alternative 
to photon cavities in many low- frequency, low-power experiments. 
The advantages of acoustic devices is the ease of operation and integ- 
ration into higher- order electronic systems, compact size, natural 
way of combination of mechanical motion and electromagnetic 
environment, as well as artificial lattice impurities. In particular, 
these systems can be used as parts of quantum hybrid systems in 
miliKelvin temperatures integrating quantum circuit electrody- 
namics and cavity elastodynamics. Moreover, unlike nano-scale 
mechanical resonators, acoustic cavities are a macroscopic objects 
of considerable effective mass varying in the range of 0.1-10 milli- 
grams for modes between 100 and 10 MHz. This, can be a promising 
experiment on unifying quantum mechanics and general relativity in 
a single experiment 7 . 

Apart from fundamental physics, cryogenic acoustic cavities can 
be used for building ultrastable frequency sources for variety of 
applications. 

Methods 

The device under test is embedded inside an oxygen-free copper block fixed on the 
end of the cold finger of a Helium pulse-tube cryocooler 2 . The cold finger is tem- 
perature-controlled, and residual fluctuations of temperature are within ±2 mK 
around the setting temperature. 

The device characterisation is performed using the impedance analyzis method: 
measurements of the displacement current through the resonator thickness while 
applying alternating potential difference between electrodes. Their ratio is the com- 
plex impedance that is measured for a range of frequencies around a resonance. This 
characteristic is used to calculate a phonon resonance bandwidth. For accurate 
measurements of intrinsic properties of a device under test (DUT), the method 
requires a calibration procedure in order to compensate any influence of a long 
connecting cables. For this purpose, three identical extra cables follow the same path 
to be ended respectively by an open-circuit, a short circuit and a 50 Ohms load. 
Measurements of these etalons helps to cancel cable influence for certain values of the 
analyzer output power, frequency span and sweep time. 

The excitation power is kept as low as possible in order to avoid nonlinear and 
thermal effects. Power dissipation inside the resonator is estimated to vary from 
0. 1 nW to 3 nW, depending on a mode coupling and quality factor. Nonlinear effects 
of the investigated devices are due to lattice anharmonicity. They strongly depend on 
quality factor of a mode and can be noticeable at powers as low as 10 nW. 

Special attention should be paid to the accuracy when characterising modes with 
extremely high quality factors due to extreme narrowness of bandwidths (below 
100 mHz). The frequency span of the analyzer is adjusted to get the maximum benefit 
of the number N of measurement points. The minimum available span of 0.5 Hz for 
401 measured points provides resolution of about 1.25 mHz. The minimum sweep 
time required to avoid all transient effects is (N — 1) X x where x is phonon relaxation 



time of a mode of frequency/ It is given as x =- 
factor, /is the resonant frequency. n J 

Based on the specifications of the measurement equipment and the post-proces- 
sing procedure, the measurement uncertainty of the intrinsic (unloaded) quality- 
factor Q is evaluated as: 



dQ 



2Q 



df 
/' 



(8) 



which demonstrates importance of the fractional frequency stability of both the 
measurement apparatus frequency source and the DUT itself. The latter stability is 
improved by locking to an ultra- stable frequency source, such as a Hydrogen maser. 
The long term (due to requirement for large sweep time) frequency stability of the 

resonator mode is determined predominantly by its temperature sensitivity and y- 

does not typically exceed 10" 10 providing measurement uncertainty below a few 
percent for all characterised modes. 
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